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Abstract 

There is a fundamental difference between the classical expression for the retarded electromag- 
netic potential and the corresponding retarded solution of the wave equation that governs the 
electromagnetic field. While the boundary contribution to the retarded solution for the potential 
can always be rendered equal to zero by means of a gauge transformation that preserves the Lorenz 
condition, the boundary contribution to the retarded solution of the wave equation governing the 
field may be neglected only if it diminishes with distance faster than the contribution of the source 
density in the far zone. In the case of a source whose distribution pattern both rotates and travels 
faster than light in vacuo, as realized in recent experiments, the boundary term in the retarded 
solution governing the field is by a factor of the order of R 1 / 2 larger than the source term of this 
solution in the limit that the distance R of the boundary from the source tends to infinity. This 
result is consistent with the prediction of the retarded potential that part of the radiation field 
generated by a rotating superluminal source decays as -R -1 / 2 , instead of i? -1 , a prediction that is 
confirmed experimentally. More importantly, it pinpoints the reason why an argument based on a 
solution of the wave equation governing the field in which the boundary term is neglected (such as 
appears in the published literature) misses the nonspherical decay of the field. 
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I. INTRODUCTION 



Scientific investigation of the electromagnetic field generated by a charged particle that 
moves faster than light began with a largely ignored article by physicist and mathematician 
Oliver Heaviside in 1887^ and is the subject of several papers by Sommerfeld in 1904 and early 
1905.^ The publication of the special theory of relativity in June 1905,^ however, discouraged 
further work as one of its tenets is, of course, that any known particle that has a charge also 
has a rest mass and so is barred from moving faster than light. Moreover, no source that 
moves faster than the wave speed can be pointlike, for this results in infinitely strong fields 
on the envelope of the emitted wave fronts. It was not until the early 1970s that Bolotovskii 
and Ginzburg^El pointed out that although special relativity precludes massive particles 
from moving superluminally, patterns of distribution of extended sources can move faster 
than light as a result of the coordinated motion of their constituent charged particles, which 
themselves remain subluminalJ 6 * 7 * 8 * 9 ^ The charge separation resulting from such coordinated 
motion gives rise to a polarization current whose distribution pattern can move faster than 
light in vacuo. The way in which this current acts as a source of radiation is immediately 

apparent from the Ampere-Maxwell equation (SI units), 

__ TT T c^D dE dP 

VXH = J + ^ = J + eo ^ + ^ (1) 
in which we see that the polarization current dP/dt contributes to the magnetic field H in 
just the same way as the current J of free charges. However, as dP/dt is not carried by 
massive particles, it is not limited to subluminal speeds. We should clarify that while the 
source distribution is superluminal, the emitted radiation (as any other) travels at the speed 
of light. 

Extended sources of electromagnetic radiation whose distribution patterns move faster 
than light in vacuo have been experimentally realized by several groups^*- 7 * 11 * 12 ! and have 
many potential applications in science and technology.^ Furthermore, they may be respon- 
sible for the extreme properties of the electromagnetic radiation received from astronomical 
objects such as pulsars (rapidly spinning, highly magnetized neutron stars) .H£l 14 l 15 l Methods 
of dealing with the asymptotic behavior of the radiation emitted by such sources, however, 
are largely neglected in standard texts.^ More dangerously, certain textbook formulae, de- 
rived in the context of stationary or subluminally moving sources, cannot be used in treating 
unusual cases such as those in which the radiation remains focused in the far zonePSDZl An 
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example is this expression for the far-zone magnetic field of a localized current with the 
density j: 

B(x P ,f P )~± /d 3 x^L. (2) 
c J |x P - x| 

Here, (xp,i P ) and (x, t) are the space-time coordinates of the observation point and the 
source points, respectively, and the square brackets denote the retarded value of V X j. (Note 
that, in contrast to the free-charge current J of Eq. (JT]), this j is a generalized current that also 
includes dT)/dt.) The derivation of Eq. ^ involves the neglect, in the far zone, of a boundary 
term containing the gradient of the radiation field. In the case of a conventional source, this 
term decays more rapidly with distance than the integral that remains in Eq. g. This is not 
the case, however, for the radiation field of a polarization current whose distribution pattern 
has an accelerated motion with a speed exceeding that of light, cJ 17 l 18 l 19 l In the next section, 
we describe such a source, the emission from which consists of a collection of narrowing 
subbeams for which the absolute value of the gradient of the radiation field B increases (as 
Rp 2 ) with the distance R-p from its source. 

Since the exact form of Eq. ^ entails an additional boundary term that depends on 



the gradient of the field in the far zone [Eq. (13), below], the increase in the magnitude 
of the gradient of the present radiation field with distance renders the boundary term in 
the retarded solution of Maxwell's equation for B dominant over its source term. The 
inadequacy of Eq. ^ in describing the radiation field generated by a superluminal source 
lies in the neglect of this boundary term. 

The remainder of this paper is organized as follows. Section II describes the 
superluminally-rotating source distribution, and a practical experimental realization of such 
a source. In order to make the later discussion more comprehensible, Section III considers 
a small volume element of this source and discusses the superposition of multiple retarded 
times that render the field of a point source divergent. Section IV generalizes this discussion 
to an extended source and describes the morphology of the resulting radiation "beam" , pro- 
viding a geometrical argument as to why the boundary term missing from Eq. 2 dominates 
in the far field, and Section V describes supporting experimental data. Sections VI and VII 
contain the main, substantive point of the current work; while the boundary contribution to 
the retarded solution for the potential can always be rendered equal to zero by means of a 
gauge transformation, the boundary contribution to the retarded solution of the wave equa- 
tion governing the field cannot be neglected in the case of a superluminal source. Finally 



Section VIII gives a brief summary. 



II. A SUPERLUMINALLY ROTATING SOURCE DISTRIBUTION 

Here we consider a particular moving source distribution for which Eq. (J2j) is not an 
adequate description of the radiation field in the far zone and its experimental realization. 
The findings of Refs. PfTT] fully clarify how this case should be treated in order to obtain 
results consistent with experimental data. In these papers, we solved the inhomogeneous 
wave equation governing the electromagnetic potential [Eq. ([7]) below] in unbounded space, 
under null initial conditions, with a polarization current density j = dP / dt for which 

P r ,ip,z(r, z, t) = s r>(PiZ (r, z) cos(m^) cos(fit), -tx < <p < n, (3) 

where 

(p = ip — tot. (4) 

Here, P r ,ip,z are the components of the polarization P in a cylindrical coordinate system based 
on the axis of rotation, u and Q are angular frequencies, s(r, z) is an arbitrary vector func- 
tion with a finite support in r > c/u, and m is a positive integer. For fixed t, the azimuthal 
dependence of the polarization ^ along each circle of radius r within the source is the same 
as that of a sinusoidal wave train, of wavelength 27rr/m, whose m cycles fit around the cir- 
cumference of the circle smoothly. As time elapses, this wave train both propagates around 
each circle with the superluminal velocity ruo and oscillates in its amplitude with frequency 
Q. This is a generic source: one can construct any distribution with a uniformly rotating pat- 
tern, P r ,tp, z (r, 0, z), by the superposition over m of terms of the form s n ^ )Z (r, z, m) cos(m(p). 
For example, the polarization-current patterns responsible for the pulsar radiation visible 
from Earth are thought to be animated by the superluminal rotation of the neutron-star 
core's intense magnetic field through the pulsar's plasma atmosphere^ one can envisage 
representing this rotating pattern by summing such terms. Moreover, Eq. (3) corresponds 
to a laboratory-based source that has been used in experimental demonstrations of some of 
the phenomena described belowP^ 

The experimental apparatus in question (see Fig. [T]) consists of a continuous strip of 
dielectric material on top of which is placed an array of metal electrodes; underneath is a 
continuous ground plate. Each upper electrode is connected to an individual amplifier, so 
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that by turning the amplifiers on and off in sequence one can apply voltages that generate a 
polarized region that moves along the dielectric with an arbitrarily high speed. In practice, 
the dielectric is a strip of alumina 10 mm thick and 50 mm across, corresponding to a 10° 
arc of a circle of average radius 10.025 m. Above the alumina strip, there are 41 upper 
electrodes of mean width 42.6 mm, with centers 44.6 mm apart 
The voltage Vj applied to the jth electrode is of the form 



where rmo and f2 are as defined in Eq. (|3]): the first cosine gives rise to the polarization- 
current wave that propagates along the dielectric and the second to a modulation of this 
wave. The source speed v is determined by the phase difference between the oscillations of 
neighboring electrodes; given the dimensions of the electrodes of the experimental machine, 
a superluminal speed v > c is obtained for At < 148.8 ps. In the experiments reported in 
Refs. [BE!, mu/(2Tr) and Q/(2tt) had the values 552.645 MHz and 46.042 MHz, respectively, 
and the source was run at v/c = 0.875, 1.064, 1.25 and 2.00 (see Refs. [6|7] for a more 
complete description). 

III. THE FIELD GENERATED BY A CONSTITUENT VOLUME ELEMENT OF 
THE SOURCE 

A superluminal source is necessarily volume-distributed P However, its field can be built 
up from the superposition of the fields of its small, moving constituent volume elements. 
Figure |2](a) shows that the waves generated by such a volume element of a rotating super- 
luminal source possess a cusped envelope and that, inside the envelope, three wave fronts 
pass through any given observation point simultaneously. This reflects the fact that the field 
inside the envelope receives simultaneous contributions from three distinct values of the re- 
tarded time [see Fig. |2{c)]. On the cusp of the envelope, where the space-time trajectory of 
the source is tangent to the past light cone of the observer [Fig. (2^d)], all three contributions 




toward the value of the field coalesce ! 8 * 21 * 22 ^ 

On this cusp (caustic), the source approaches the observer with the speed of light and 
zero acceleration at the retarded time, i.e. dR(t)/dt = —c and d 2 R(t)/dt 2 = 0, where 
R(t) = |x(t) — xp| is the distance between the source point x(i) and the observation point 



Vj = Vq cos[mu;(j At — t)) cos(fit) 



(5) 
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FIG. 1: Left: experimental animation of a superluminal polarization current, (a) A simplified 
dielectric solid containing negative (0) and positive (©) ions. In (b), a spatially- varying electric 
field has been applied, causing the positive and negative ions to move in opposite directions; a 
finite polarization P has therefore been induced. If the spatially-varying field is made to move 
along the direction of the arrow, the polarized region moves with it. (c) Schematic side view 
of a practical superluminal source, showing metal electrodes above a strip of dielectric (shaded 
region) and a ground plate below it. "0" indicates that there is no voltage on that particular upper 
electrode; the symbol + indicates a positive voltage applied to the upper electrode. The voltage 
on the electrodes produces a finite polarization of the dielectric (darker shading) . (d) By switching 
the voltages on the electrodes on and off, the polarized region (darker shading) can be made to 
move along the dielectric, (e) Top view, showing the curvature of the dielectric (lighter shaded 
region). The curvature introduces centripetal acceleration in the moving polarized region. Note 
that the electrodes (black shading) cover only part of the top surface of the dielectric. Right: a 
photograph of the experimental apparatus showing a curved strip of dielectric material (alumina) 
sandwiched between a copper ground plate (below) and 41 electrodes, each connected to its own 
amplifier (left). 

xp. As a result, the interval of emission time for the signal carried by the cusp is much 
longer than the interval of its reception timeP^ Fig. |2^d) represents this diagrammatically; 
rather than merely intersecting at one or three points, the space-time trajectory of the source 
and the observer's light cone become tangent to one another. This essentially instantaneous 
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(b) (c) (d) 



FIG. 2: (a) Envelope of the spherical wave fronts emanating from a superluminally moving source 
element (S) in uniform circular motion. The lighter circles are Huygen's wavelets emitted by the 
source as it traverses its circular orbit, designated by the larger of the two dotted circles. The 
heavy curves show the cross section of the envelope of the wavelets within the plane of the orbit of 
the source. The smaller of the dotted circles represents the light cylinder r = c/lo. (b), (c) and (d) 
Space-time {i.e. ct versus distance x) diagrams showing the intersection of the trajectory of the 
source point S with the past light cone of the observation point P when P lies outside (b), inside 
(c), and on the cusp of (d) the envelope of wave fronts. 

reception of contributions from an extended period of emission time represents focusing 
of the radiation in the time domain; it has been described as temporal focusing and is 
the subject of two experimental papersP^ It is this temporal focusing that leads to the 
divergence of the field of a point source on the cusp of the envelope of wave fronts and is 
ultimately responsible for the component of the radiation from an extended source whose 
intensity decays as 1/distance, rather than the inverse square law that applies to all other 
sources^^, consequently, this part of the radiation will dominate in observations made from 
large distances. 

A three-dimensional view of the envelope of wave fronts and its cusp is shown in Fig. |3j 
The two sheets of the envelope, and the cusp along which these two sheets meet tangen- 
tially, spiral outward into the far zone. In the far zone, the cusp lies on the double cone 
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Light cylinder 



(a) (b) 
FIG. 3: Three-dimensional views of the envelope (a) and its cusp (b). 

9 P = arcsin[c/(ra>)], 9 P = n — axcsinfc/ (ru;)], where (Rp,0p,<pp) denote the spherical polar 
coordinates of the observation point P. Thus, a stationary observer in the polar interval 
arcsin[c/(rc<j)] < Op < it — arcsin[c/(ra;)] receives recurring pulses from each volume element 
as the envelope rotates past him/herP 

Figure [4] shows the radiation field generated by the rotating source element S on a cone 
close to the cusp, just outside the envelope. Not only does the spiraling cusp embody a 
recurring pulse, but the plane of polarization of the radiation swings across the pulsep^ as 
does the radio emission received from pulsars 



IV. THE FIELD GENERATED BY THE ENTIRE VOLUME OF THE SOURCE 

The dominant contribution towards the field of an extended source comes from a thin 
filamentary part of the source that approaches the observer, along the radiation direction, 
with the speed of light and zero acceleration at the retarded timePFor an observation point P 
in the far zone with the coordinates (Rp, 6p, fp), this filament is located at r = (c/oj) csc#p, 
ip = Lfp + 3tt/2 and is essentially parallel to the rotation axis (Fig. [5]). The collection 
of cusps of the envelopes of wave fronts that emanate from various volume elements of 
the contributing filament form a subbeam whose polar width is nondiffracting: the linear 
dimension of this bundle of cusps in the direction parallel to the rotation axis remains the 
same at all distances from the source, so that the polar angle 59p subtended by the subbeam 
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-5 0, 5 

x p co/c 

FIG. 4: Polarization position angles and field strengths on the cone Op = n/12 outside the envelope 
of wave fronts for a source element with ruj = 2. Here, the field vectors are projected onto the 
plane (ccp,yp) of the source's orbit. 

decreases as -Rp -1 with increasing Rp (Fig. [5] and Ref. [9]). 

In that it consists of caustics and so is constantly dispersed and reconstructed out of 
other waves, the subbeam in question radically differs from a conventional radiation beam 
(see Appendix D of Ref. The narrowing of its polar width (as -Rp -1 ) is accompanied by 

a more slowly diminishing intensity (an intensity that diminishes as Rp' 1 instead of Rp~ 2 
with distance), so that the flux of energy across its cross sectional area remains the same 
for all Rp.- As we will see in the next section, this slower rate of decay of the emission from 
a superluminally rotating source has been confirmed experimentally upto several hundred 
Fresnel distances (Rayleigh ranges) 

The contributing part of an extended source (the filament that approaches the observation 
point with the speed of light and zero acceleration) changes as the source rotates (see Fig. [5]). 
In the case of a turbulent plasma with a superluminally rotating macroscopic distribution, 
therefore, the overall beam within which the narrow, nonspherically decaying radiation is 
detectable would consist of an incoherent superposition of coherent, nondiffr acting subbeams 
with widely differing amplitudes and phases (similar to the train of giant pulses received 
from the Crab pulsar).^ 
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FIG. 5: Schematic illustration of the light cylinder r = c/uj, the filamentary part of the source 
that approaches the observation point with the speed of light and zero acceleration at the retarded 
time, the orbit r = c/(wsin#p) of this filamentary source, and the subbeam formed by the bundle 
of cusps that emanate from the constituent volume elements of this filament. 

The overall beam occupies a solid angle whose polar and azimuthal extents are inde- 
pendent of the distance Rp. It is detectable within the polar interval arccos^r^/c) -1 ] < 
\6p — ?r/2| < arccos[(r u to , /c)~ 1 ], where [n,r u ] denotes the radial extent of the superluminal 
part of the source. The azimuthal profile of this overall beam reflects the distribution of the 
source density around the cylinder r = c/(u sin dp), from which the dominant contribution 
to the radiation arisesP 

Because the subbeams that constitute the overall beam are narrower the farther away 
from the source they are detected, the absolute value of the gradient of the radiation field 
associated with them increases with distance: it has been explicitly shown in Ref. [17] 
that the gradient of the nonspherically decaying field generated by the superluminal source 
described in Eq. ^ has a magnitude that increases as R P 7/2 with increasing; Rp. 

V. EXPERIMENTAL AND OBSERVATIONAL CONFIRMATION THAT THE 
NONSPHERIC ALLY-DECAYING TERM DOMINATES IN THE FAR FIELD 

The issue under discussion in this paper is the dominance of the boundary term neglected 
in Eq. (2) in the far field. Experimental demonstrations of the far-field dominance of the 
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FIG. 6: (a) and (b): Ratio of the detected intensity with the experimental source (Fig. [T| running 
superluminally (at v/c = 1.064) to that with the source running subluminally (v/c = 0.875) versus 
distance Rp. For both data sets the detector was moved along a path close to the direction of 
the cusp expected for v/c = 1.064. Data are points; the line is a straight line through the origin, 
which implies that the intensity in the superluminal case would decline in free space as l/i?p (after 



boundary term focus on detecting emissions whose intensity decays with distance as l/Rp, 
rather that the more conventional inverse-square law. 

Fig. [6] shows a ground-based demonstration^ using the machine depicted in Fig. [I] The 
experiment consisted of moving a detector along the expected cusp direction for a source 
speed of v/c = 1.064 and recording the intensity as a function of distance. The source speed 
was then set to a subluminal velocity (i.e. no cusp formation- intensity follows the inverse- 
square law with distance) and the experiment repeated with the detector following the same 
path. The ratio of the detected intensity for v/c = 1.064 to that for the machine running at 
subluminal speed increases linearly with distance, showing that the cusp radiation intensity 
is decaying more slowly with increasing distance, in line with expectationsP 

For demonstrations covering much larger distances, one must turn to astronomical ob- 
servations. In a recent paper,®! ^ ne Maximum Likelihood Method 54 ' has been applied to 



Ref. (6)). 
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FIG. 7: Demonstration that the far-field emission of pulsars is dominated by a component that 
decays nonspherically with distance (after Ref. (33)). (a) The relative convergence error from the 
Maximum Likelihood Method applied to 971 intragalactic pulsars (from the Parkes Multibeam 
Survey) as a function of the trial exponent retrial- There is a clear minimum at retrial = 1 5 strongly 
suggesting that the flux at 1400 MHz, S'1400 oc l/Rp. (b) The inferred luminosity function with 
raw pulsar data (taken from the Parkes Multibeam Survey) plotted as large circles, and the fitted 
luminosity function as small diamonds. It is clear that there is a very close correspondence between 
fit and data. 

observational data from 971 pulsars to deduce the likely power law with which their flux 
S (a quantity proportional to intensity) falls off with distance. The Maximum Likelihood 
Method determines the probable luminosity function (i.e. the distribution of pulsars with 
respect to their radiated power) whilst making allowances for missing observations due to 
the low sensitivity of Earth-bound instruments.^ The statistical success of the method is 
measured using a relative convergence error; to find the exponent n in the relationship 
S oc l/Rp, the calculation is repeated for several different values of a trial exponent, n tr iai 
until the value that gives the smallest relative convergence error is located. As can be seen 
in Fig. [7j the pulsar data are most convincingly fitted by a flux S oc 

In summary, experimental and observational data suggest very strongly that the boundary 
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term neglected in Eq. (2) dominates in the far field, resulting in a component of the radiation 
that is nonspherically decaying, with an intensity proportional to 1/Rp. 



VI. FUNDAMENTAL ROLE OF THE RETARDED POTENTIAL 

While the results of an analysis based on the retarded solution to the wave equation 
governing the field depend crucially on the boundary conditions satisfied by the field at 
infinity,^ if we instead base our analysis on the retarded potential, we require no corre- 
sponding explicit knowledge of the value of the potential in the radiation zone. By using 
the field derived from the retarded potential [Eq. (JlT|], we will evaluate the boundary term 



in the retarded solution to the wave equation governing the field [Eq. (13)] and show that, 
far from being equal to zero, as assumed in Refs. |16f25f26|l27f28f29| . this boundary term 
constitutes the dominant contribution to the value of the field in the far zone. There is 
no discrepancy between the results obtained from the retarded solution for the potential 
and the retarded solution for the field once the boundary term in the solution to the wave 
equation governing the field is retained. Furthermore, this analysis is robust with respect to 
choice of integration boundaries (see the Appendix). 

We first describe how the boundary contribution to the retarded solution for the potential 
can always be made equal to zero, irrespective of the source motion. 

In the Lorenz gauge, [i.e., the choice of a set of potentials (A, A ) that satisfy the Lorenz 
condition V-A + c~ 2 dA°/dt = 0], the electromagnetic fields 

E = -VpA° -dA/d(ct P ), B = V P xA, (6) 

are given by a four-potential A 11 that satisfies the wave equation 

1 f) 2 A^ Air 

v 2 ^-4^r = ~-^ A* = 0,.-,3 ) (7) 
cr at 1 c 

where A / c and j° / c are the electric potential and charge density and A^ and for // = 1, 2, 3 
are the Cartesian components of the magnetic potential A and the current density j.^ The 
solution to the initial-boundary value problem for Eq. ^ is given by 

A"(x P , t P ) =- f V dt [ d 3 x fG + — r dt [ dS • (GVA» - A^VG) 

C Jo Jv 47r Jo Jy, (o\ 

-4^| d H # ¥- G ¥),o' 
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in which G is the Green's function and £ is the surface enclosing the volume V (see, e.g., 
page 893 of Ref. [30]). 

The potential arising from a general time-dependent localized source in unbounded space 
decays as Rp 1 when Rp = |xp| — > oo, so that for an arbitrary free-space potential the second 
term in Eq. (fsl) would be of the same order of magnitude (~ -Rp 1 ) as the first term in the 
limit that the boundary X tends to infinity. However, even those potentials that satisfy the 
Lorenz condition are arbitrary to within a solution of the homogeneous wave equation: the 
gauge transformation 

A^A + VA, A°^A°-dA/dt (9) 

preserves the Lorenz condition if V 2 A — c~ 2 d 2 A/dt 2 = 0.^ One can always use this gauge 
freedom in the choice of the potential to render the boundary contribution (the second term) 
in Eq. ^ equal to zero, since this term, too, satisfies the homogeneous wave equation. Under 
the null initial conditions A^\ t=0 = (<9v4 M / 'dt) t =o = assumed in this paper, the contribution 
from the third term in Eq. ^ is identically zero. 

In the absence of boundaries, the retarded Green's function has the form 

G(^ P ,t P ) = S{tP ~ t ~ R/C \ (10) 

where 5 is the Dirac delta function and R is the magnitude of the separation R = xp — x 
between the observation point xp and the source point x. Irrespective of whether the 
radiation decays spherically (as in the case of a conventional source) or nonspherically (as 
applies for a rotating superluminal source — see Section |h))J^^ therefore, the potential A^ 1 
due to a localized source distribution that is switched on at t = in an unbounded space, 
can be calculated from the first term in Eq. pi): 



A"(x P ,tp) = c- 1 J ^x&tf{^t)5{t v ~t-R/c)/R, (11) 

i.e., from the classical expression for the retarded potential. Whatever the Green's function 
for the problem may be, in the presence of boundaries, it approaches the expression in 



Eq. (10) in the limit where the boundaries tend to infinity. Hence one can also use this 



potential to calculate the field on a boundary that lies at large distances from the source. 
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VII. RETARDED SOLUTION OF THE EQUATION GOVERNING THE FIELD 



We now return to the case of the field and show that an analogous assumption about the 
boundary contribution may not be made. Consider the wave equation 

1 <9 2 B 4tt „ . 



governing the magnetic field; Equation ( 12 ) may be obtained by simply taking the curl of 



the wave equation for the vector potential [Eq. ([7| for n = 1, 2, 3]. We write the solution to 



the initial-boundary value problem for Eq. (12), in analogy with Eq. <\8h, as 



J B fc (x P ,tp)=- r dt [ d 3 x{Vxj) k G + ^- r dt [ dS • {GVB k - B k VG) 
c Jo Jv 47r Jo Jt, 

1 ' A fn«.c»4 



(13) 



Arte 2 J v V dt dt / t=o 
where = 1,2,3 designate the Cartesian components of B and Vxj. 

Here, we no longer have the freedom, offered by a gauge transformation in the case 
of Eq. ([8]), to make the boundary term zero, nor does this term always decay faster 
than the source term, so that it could be neglected for a boundary that tends to infin- 
ity, as is commonly assumed in textbooks (e.g., page 246 of Ref. [U5]) and the published 
literature. 1 25 1 26 1 27 1 2 ^ 1221 The boundary contribution to the retarded solution of the wave equa- 



tion governing the field [the second term on the right-hand side of Eq. (13)] entails a surface 
integral over the boundary values of both the field and its gradient. For the rotating su- 
perluminal source (3), where the gradient of the field increases as Rp 2 over a solid angle 
that decreases as -Rp 4 , this boundary contribution is proportional to Rp 1 ^ 2 (see Ref. [PT]). 
Not only is this not negligible relative to the contribution from the source term [the first 



term on the right-hand side of Eq. (13)], but the boundary term constitutes the dominant 
contribution to the radiation field in this case^ 

If one ignores the boundary term in the retarded solution of the wave equation governing 
the field, as in Refs. [16,25.26,27,28.29] , one obtains a different result, in the superluminal 
regime, from that obtained by calculating the field via the retarded potential) 9 * 31 * 32 ! This 
apparent contradiction stems solely from having ignored a term in the solution to the wave 

1/2 

equation that is, by a factor of the order of Rp , larger than the term that is normally kept. 
The contradiction disappears once we take the neglected term into account: the solutions 
to the wave equations governing both the potential and the field predict that the field of a 
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rotating superluminal source decays as R P as Rp tends to infinity, a result that has been 
demonstrated experimentally.^ 

We note, furthermore, that the representation 

A(x P ,tp) = - [d*x^\, (14) 
cj |x-x P | 

of the retarded potential is differentiable as a classical (as opposed to generalized) function 
only in the case of a moving source whose speed does not exceed that of the waves it 
generates. The steps, familiar from the subluminal regime, one takes to derive Eq. (J2| 



by differentiating Eq. (14) as a classical function are not mathematically permissible when 
the moving source has volume elements that approach the observer with the wave speed 
and zero acceleration at the retarded timeP^ Contrary to the usual assumption] 25 * 26 * 27 * 28 * 29 ! 
the retarded distribution of the density of a moving source is not necessarily smooth and 
differentiable if its rest-frame distribution is. The retarded distribution of a rotating source 
with a moderate superluminal speed is in general spread over three disjoint volumes (differing 
in shape from each other and from the volume occupied by the source in its rest frame) whose 
boundaries depend on the spacetime position of the observer. The limits of integration in 



Eq. (14), which delineate the boundaries of the retarded distribution of a localized source, 
are not differentiable functions of the coordinates of the observer at those points on the 
source boundary that approach the observer, along the radiation direction, with the speed 
of light at the retarded timeP^*^ In the superluminal regime, derivatives of the integral 
representing the retarded potential are well-defined only as generalized functionsP 

VIII. SUMMARY AND CONCLUSIONS 

We must stress that we have not calculated the radiation field directly from the solution 
to the wave equation governing the field (the solution whose boundary term is normally 
neglected) j 25 l 26 l 27 l 28 l 29 l We have instead first solved the wave equation governing the potential 
(whose solution has no boundary term) and then used this solution to evaluate the neglected 



term in the exact version of the retarded solution for the field [Eq. (13)]. What one obtains 
by including the boundary term in the retarded solution to the wave equation governing the 
field is merely a mathematical identity; it is not a solution that could be used to calculate the 
field arising from a given source distribution in free space. Unless its boundary term happens 
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to be small enough relative to its source term to be neglected, a condition that cannot be 
known a priori, the solution in question would require that one prescribe the field in the 
radiation zone (i.e., what one is seeking) as a boundary condition. Thus, the role of the 
classical expression for the retarded potential in radiation theory is much more fundamental 
than that of the corresponding retarded solution of the wave equation governing the field. 
The only way to calculate the free-space radiation field of an accelerated superluminal source 
is to calculate the retarded potential and differentiate the resulting expression to find the 
fido™ (see also Refs. |5TI32| ). 
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APPENDIX A: A NOTE ON INTEGRATION BOUNDARIES 

It has been suggested^ that our demonstration of the nonspherical decay of the field of 
a rotating superluminal source can be dismissed on the grounds that the boundary term in 



Eq. (13) must be identically zero. However, the claim that the boundary contribution to 
the retarded solution of the wave equation governing the field is exactly zero even when the 
boundary encloses a source violates the very foundations of diffraction theory: For a volume 



V in which no sources are present, Eq. (13) reduces to 

£? fc (x P ,tp) = i- / P dt I dS • (GVB k - B k VG), (Al) 

J J s 

under the null initial conditions B k \ t= Q = (dB k /dt) t=0 = 0. As in the customary geometry for 
diffraction, the closed surface E can consist of two disjoint closed surfaces, S inner and E outer 
(e.g., two concentric spheres), both of which enclose the source (see Fig. 10.7 of Ref. [To]). 
If the observation point does not lie in the region between Ei nner and E oute r, i-e., lies outside 



the closed surface E, then the composite surface integral in Eq. (|A1|) vanishes: 

fc VG) 

(A2) 



P dt / dS ■ (GVB k - B k VG) = [ ? dt ( / + / 1 dS • (GVB k - B k VG) 



= 



outer 



(see Ref. [30]). Under no circumstances, however, would the integrals over E inner or E outer 
vanish individually if these surfaces enclose a source, i.e., if there is a nonzero field inside 
E^ner; n or does the invariance of the values of these integrals under deformations of E in ner 



and S outor have any bearing on whether they are nonvanishing or not™ Equation (Al) 
forms the basis of diffraction theory If the surface integrals over E imi er and E outer were to 
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vanish individually, as claimed in Ref . [29J , the diffraction of electromagnetic waves through 
apertures on a surface enclosing a source would be impossible (see Sec. 10.5 of Ref. [TB]). 
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